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Systcas  are  Modeled  in  order  to  understand  and  explain  thea  better  and  as  a 
prelude  to  action.  Aircraft  dynaaics,  for  exaaple,  may  identified  so  that 
better  designs  can  be  aade,  or  so  that  adaptive  control  actions  can  be  taken. 

Our  attention  in  this  continuing  study  has  been  directed  at  aspects  of  estiaati  >n 
and  identification  that  are  connected  with  systea  understanding.  This  study  ha*^ 
been  aiacd  at  state  estiaation  and  paraaeter  identification  for  a  new  class  of 
Models,  causal  functional  equations,  which  describe  wave  propagation  in  layered 
aedia  systcoM.  These  oxKlels  are  applicable  to  diverse  areas,  such  as  reflectio 
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I .  INTRODUCTION 


SysCema  are  nodeled  in  order  to  underaCand  and  explain  thea  better 
and  aa  a  prelude  to  action.  Aircraft  dynaaica,  for  exaaple,  may  be  iden¬ 
tified  80  that  better  designs  can  be  Mde,  or  so  that  adaptive  control 
actions  can  be  taken.  Our  attention  in  this  continuing  study  has  beer 
directed  at  aspects  cf  estiouitior  and  identification  that  art  connected 
vltb  systee  understanding.  This  study  has  been  aimed  at  state  estimation 
and  parameter  identification  for  a  new  class  of  models,  causal  functional 
equations,  which  describe  wave  propagation  in  layered  media  systems.  These 
models  are  applicable  to  diverse  areas,  such  as  reflection  seisonlogy, 
transmission  lines,  speech  processing,  optical  thin  coatings  and  EM  problems. 


II.  RESEARCH  PROGRESS 

1.  An  imiortant  special  case  rf  a  causal  functional  eqcatlon  fCFE) ,  occurs 
when  all  onm-way  travel  times  are  equal.  In  this  case  the  uniform  CFE  (UCFE) 

is 

x{t  +  t)  -  A  x(t)  +  b  m(t)  (1) 

with  initial  values 

x(o)  o  c  [0,  t)  A  J"  (2) 

Recognising  that  any  time  t  (t  (  R)  can  be  expressed  as 

t  ■  t*  +  Mt  %rttere  t*(  J  and  M  is  an  integer,  (3) 

we  have  shown  [1]  that  the  solution  to  UCFE  (1)  is 

k 

X  [t*  +  (k+l)T]  -  x(t’>  +  m(t'+  It)  (4) 

1-0 


where  t'  c  J^,  k  -  0,  1,  2,  ...,  and  t  -  t*  +  (k+l)T. 

Equation  f4>  explicitly  shows  how  the  state  at  any  time  t  -  t’  +  (k+l)T 
depends  or  an  initial  condition  x(t’)  and  the  input  m.  It  is  of  interert  to 
note  that  x{t)  depends  only  on  a  single  element  of  the  initial  values  x(o) 

(o  (  J)t  namely  and  a  finite  number  of  point  values  of  m.  This  shows 

that  the  solution  totheuniform  causal  functional  state  equation,  although 
contlnuous-tlaw  in  nature,  derives  its  values  in  a  discrete-time  fashion  for  a 
given  fixed  value  of  t*  f  J.  Of  course,  there  are  an  vncountable  number  of 
points  In  J";  hence  we  can  imagine  x(t)  as  being  generated  by  an  uncountable 
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When  we  sinulaCe  our  results  on  a  digital  coaputer,  conputatlons  are 
made  every  T  sec.  at  discrete  tine  points.  Consequently,  on  a  digital 
conputer,  x(t)  Is  generated  by  a  finite  ruaber  of  dlscretc-tlme  rysteas 
vhlch  o(>erare  In  parallel. 

2.  We  have  derived  [ij  the  alnlmua-varlance  state  estimator  for  UCFE 

x(t  +  t)  -  A  x(t)  +  B  m(t)  +  w(t)  (5) 

and  Its  associated  measureaent  equation 

i(t)  -  H  x(t)  -h  n(t)  (6) 

Let  x(t)  denote  the  min Inum-var lance  estimate  of  ]((t)  based  on  all 
measureaents  In  {^(1):  0^Xj<t,  tcR}.  We  have  shown  that  for  any  fixed 

t*  c  J"- ■  [0,t),  ^(t),  where  t  ■  t '  +  Mx  (M  ■  1,  2,  ...),  Is  given  by  the  usual 
discrete-time  Kalman  filter  equations  (Ref.  A,  for  example)  with  t*  considered 
the  Initial  starting  time.  Of  course,  to  obtain  ^(t)  for  all  t  c  R  we  would 
need  an  uncountable  number  of  discrete-time  Kalman  filters;  but,  imposing  a 
mesh  on  J'  (with  grid  size  equal  to  data  sampling  rate)  leads  to  a  finite 
number  of  Kalman  filters  which  operate  In  parallel.  To  the  best  knowledge  of 
the  author  this  Is  the  first  estimation  theory  result  that  has  led  to  a  natural 
fora  of  parallel  dr.ta  processing. 

3.  We  have  developed  an  extended  minimum* variance  estimator  for  a Imul taneoi'.s 
estimation  of  states  and  parameters  (i.e.,  reflection  coefficients)  in  a  UCFE 
[23.  Simulation  results  were  very  disappointing.  Reasons  for  the  disappointing 
results  are  explained  in  [2]. 

4.  A  slaple  Inverse  filter  has  been  developed  [4]  to  suppress  multiple 
reflections  from  a  normal  Incidence  synthetic  seismogram.  The  filter  was 
developed  by  means  of  Mendel's  Brenner  Series  Decomposition  [3]  and  the  operator 
description  of  state  spare  mode)  of  layered  media  (Ref.  B) ,  and  Is  given  In 
terms  of  z-transfonss  as 
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In  this  equation  Y(x)  is  the  synthetic  selstsogran  measured  by  a  sensor 
located  on  the  surface.  That  sensor  receives  reflected  signals  from  a 
layered  media  due  to  the  normal  Incident  Input  M(z)  applied  at  the  same 
surface.  Parameter  r^  is  the  reflection  coefficient  of  the  surface.  ?^(z) 

Is  the  output  of  the  filter;  It  consists  of  the  primary  reflection  portion 
of  the  seismogram  and  some  residual  terms;  l.e.,  Y^(s)  •  Y^(c)  -f  y(z)*  In 

aw 

general,  the  residual  terms  y(z)  Is  relatively  quite  small,  and  Yj^(z)  Is  a 
good  approximation  of  Y^(z).  This  filter  is  especially  effective  when  r^ 

Is  relatively  large  (as  In  most  geophysical  situations)  in  which  case  y(z) 

Is  almost  negligible  compared  with  Y^(z).  The  filter  requires  knowledge  of 
the  Input  waveform  M(z) ,  surface  reflection  coefficient,  r^,  and  measured 
seismogram  Y(z).  Observe  that  (7)  represents  e  nonlinear  processing  of  the 
seismogram  Y(z). 

5.  We  have  extended  our  normal  Incidence  state  space  model  to  the  non-normal 
Incidence  case  [5].  The  non-normal  Incidence  (NNI)  state  space  model  Is 
structurally  the  same  as  the  normal  Incidence  state  space  model  except  that  It 
has  twice  as  many  state  variables.  Because  of  mode  conversion  In  non-normal 
Incidence,  the  scalar  upgolng  and  downgoing  waves  and  travel  times  In  each 
layer  as  well  as  reflection  and  transmission  coefficients  in  each  interface  are 
replaced  by  a  vector  of  upgolng  and  downgoing  waves,  a  vector  of  travel  times, 
and  matrices  of  reflection  and  transmission  coefficients. 

With  this  NNI  model,  we  are  able  to  generate  synthetic  seismograms  for  a 
plane  wave  source,  and  more  importantly,  for  a  two-dimensional  point  source. 

6.  We  have  developed  a  maximum-likelihood  procedure  for  estimating  both  thc^ 

reflection  coefficients  and  one-way  travel  times  [6]  for  a  lossless  layered 
media  system  in  which  the  layers  are  non-equally  spaced  (In  time).  It  uses  a 
state  space  aiodel  as  its  starting  point,  one  that  is  more  general  than  (1) 
since  now  t  is  different  for  each  layer  (l.e.,  ^  ''^2  ^  ^3  ^  ^k  ^ 

The  only  source  of  uncertainty  Is  measurement  noise,  n(t).  Haxlmum-llkel ihood 
estimates  of  the  par.-uaeters  are  obtained  In  a  layer-recursive  format.  In 
essence,  first  r^  and  are  determined  and  layer  1  Is  stripped  away;  then  r^ 
and  Tj  are  determined  and  layer  2  is  stripped  away;  etc.  To  the  best  of  our 


knowledge,  thie  Is  the  first  tine  that  both  reflection  coefficients  and 
travel  tines  have  been  slnultaneously  estlnated  in  an  optinal  manner. 


PUBLICATIONS  UNDER  GRANT  75-2797 


[Ij  H.  D.  Washburn  and  J.  M.  Mendel,  "Minlnun-Varlance  State  Estimation 
for  Unlforn  Causal  Functional  Equations,"  presented  at  17th  IEEE 
Conference  on  Decision  and  Control,  San  Diego,  California,  January,  1979. 

[2j  J.  M.  Mendel,  "Identification  of  Reflection  Coefficients  fron  Noisy  Data 
by  Means  of  Extended  Mlnimun  Variance  Estimators:  A  Critical  Exanlnatlon, 
presented  at  17th  IEEE  Conference  on  Decision  and  Control,  San  Diego, 
California,  January,  1979. 

J.  H.  Mendel,  "Brenner  Series  Deconpositlon  of  Solutions  to  the  Lossless 
Wave  Equation  In  Layercxl  Media,"  IEEE  Trans,  on  Geoscience  Electronics, 
Special  Issue  on  Seismic  Data  Processing,  April,  1978. 

[4]  J.  s.  Lee  and  J.  M.  Mendel,  "Suppression  of  Multiple  Reflections," 
presented  at  48th  Annual  International  Meeting  of  the  Society  of 
Exploration  Geophysicists,  San  Francisco,  November,  1978. 

[53  F.  Amlnzadeh  and  J.  M.  Mendel,  "Non-Normal  Incidence  State  Space  Model," 
presented  at  48th  Annual  International  Meeting  of  the  Society  of 
Exploration  Geophysicists,  San  Francisco,  November,  1978. 

[6]  F.  H.  Ashrafl,  "Estimation  of  Parameters  in  Lossless  Layered  Media 
Systems,"  Ph.D.  dissertation.  University  of  Southern  California, 

October, 1978. 


REFERENCES 


A.  J.  S.  Medltch,  "Stochastic  Optimal  Linear  Estimation  and  Control," 
McGraw-Hill,  New  York,  1969. 

B.  J.  M.  Mendel,  N.  E.  Nahl,  L.  M.  Silverman,  and  H.  D.  Washburn,  "State 
Space  Models  of  Lossless  Layered  Media,"  presented  at  1977  Joint 
Automatic  Control  Conference,  San  Francisco,  California,  June,  1977. 


